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OBSERVATIONS  ON  APPROXIMATE  AND  EXACT 
TREATMENTS  OF  FRACTURE  STATISTICS 
FOR  POLYAXIAL  STRESS  STATES 


In  his  original  treatment  of  the  statistics  of  fracture,  Weibull  [ l] 
introduced  a two- parameter  functional  form  for  the  relationship  between  sim- 
ple tension  and  probability  of  fracture.  He  also  showed  how  to  compute  the 
corresponding  statistics  of  failure  for  bending,  torsion,  and  other  stress 
states  that  involve  only  uniaxial  tensile  stresses.  These  applications  have 
become  well  known  and  are  widely  used  [2-4], 

In  problems  that  involve  biaxial  or  triaxial  tensile  stresses,  the  situa- 
tion is  more  complicated.  Weibull  gave,  without  formal  proof,  a procedure 
for  treating  such  problems  and  showed  how  to  apply  it  in  several  simple  cases. 
Some  investigators  have  expressed  doubts  concerning  the  rigor  of  Weibull's 
treatment  of  polyaxial  stress  states  [5,6],  and  there  are  indications  that  these 
doubts  were  later  shared  by  Weibull  himself  [7,8].  Since,  in  addition,  rather 
tedious  calculations  are  required  for  each  polyaxial  stress  state,  in  practical 
structures  that  involve  continuously  varying  stress  states,  there  is  a natural 
tendency  to  use  approximations.  One  simple  approximation,  which  constitutes, 
in  fact,  the  only  technique  suggested  for  handling  polyaxial  stress  statistics  in 
a recent  treatise  on  fracture  [9],  is  to  assume  that 

Ps  (al’  a2’  a3}  = Ps  (al}  Ps  (a2>  Ps  (a3)  (1) 

where  Pg  is  the  probability  of  survival,  and  (jj,  o^,  a 3 are  the  principal 
stresses.  This  method  was  mentioned  earlier  by  Barnett  et  al.  [5].  For 
brevity,  it  is  referred  to  here  as  the  Barnett-Freudenthal  (BF)  approximation. 

One  objective  of  the  present  investigation  is  to  investigate  the  limitations 
of  this  approach  and  thus  assist  potential  users  in  understanding  the  general 
nature  and  approximate  magnitude  of  the  errors  involved.  Another  objective 
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is  to  point  out  an  alternative  approximation  that  is  preferable  in  some 
applications.  For  the  sake  of  simplicity  and  brevity,  interest  is  focused 
primarily  on  biaxial  tension  applied  to  isotropic  materials  whose  fracture 
statistics  in  simple  tension  can  be  described  with  satisfactory  accuracy  by 
Weibull's  two-parameter  formulation: 

Pf  = 1 - exp  [-kam]  (2) 

or 

in  Pg  = in  (1  - Pf)  ■■  -kom  (3) 

where  Pj.  is  the  probability  of  fracture.  An  exact  treatment  of  the  Weibull 
theory  solution  is  also  given. 

The  BF  approximation  would  be  strictly  correct  if  all  the  cracks  were 
oriented  with  their  planes  normal  to  any  one  of  the  principal  stresses.  How- 
ever, this  is  rarely  the  case  and  certainly  does  not  occur  in  isotropic  ma- 
terials. Generally,  there  will  be  some  cracks  inclined  to  the  principal  axes 
that  will  be  fractured  by  combined  stresses  even  though  capable  of  surviving 
any  of  the  principal  stresses  applied  individually.  This  phenomenon  tends  to 
make  eq.  (1)  unconservative.  However,  there  are  also  cases  in  which  a given 
crack  would  be  fractured  by  more  than  one  of  the  principal  stresses  even  acting 
alone,  and  this  tends  to  make  the  assumption  of  noninteraction,  i.  e.  , eq.  (1), 
conservative.  We  show  how  these  opposite  trends  interact,  and  that,  gener- 
ally, the  unconservative  tendency  dominates. 

In  the  case  of  equibiaxial  tension,  eq.  (1)  reduces  to 

P_  (cr»  a)  = [P  (a,  o)]2  (4) 

S b 
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or 


in  Pg  (a,  a) 
in  Pg  (a,  o)  = 2 


(5) 


Equation  (5)  implies  that,  for  the  small  values  of  P^.  of  primary  interest 

in  Pg  (a,  a)  in  [ 1 - Pf  (a,  a)]  Pf  {a,  a) 

in  Pg  (a,  o)  = in  1 1 - P^  (a,  o)]  ~ Pf  (a,  o)  = 2 


i.  e.  , the  probability  of  failure  in  equibiaxial  tension  is  twice  that  for 
uniaxial  tension. 

As  noted  earlier,  one  would  expect,  from  general  principles,  that 
this  would  sometimes  be  an  underestimate  of  the  probability  of  equibiaxial 
fracture.  In  order  to  determine  when  and  how  much,  we  must  turn  to 
theory.  For  this  purpose,  we  select  a recent  reformation  of  weakest  link 
theory  in  which  the  flaws  are  identified  as  cracks  [ 1 0 J , and  the  polyaxial 
stress  statistics  are  derived  in  a straightforward  manner  from  the  basic 
assumptions.  It  was  noted  earlier  [8]  that,  according  to  this  theory,  when 
the  probability  of  fracture  in  simple  tension  obeys  eq.  (2),  the  probability 
for  equibiaxial  tension  becomes 


Pj.  (a,  a)  = 1 - exp 
For  integral  values  of  m,  eq. 

Pj.  (a,  a)  = 1 - exp 


T i~>  j.i\  i m r(m)  T (1.5)1 
j^-(2m  + l)mka  r(m  + 1.5)J 

(7)  becomes 

T k (m)  \ qm  1 

" (m  - 0.  5)  (m  - 1.  5)  --  0.5 


(7) 


(8) 
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Combining  eqs.  (2)  and  (8),  we  obtain 


In  P 

TTP 


s 


s 


(a>  CT)  m • 

(a,  o)  " (m  - 0.5)  (m  - 1.5)  ---  0.5 


(9) 


In  Fig.  1,  this  result  is  compared  with  the  previous  one  given  in 
eq.  (5),  and  we  conclude  that  the  BF  approximation  leads  to  the  correct 
result  for  equibiaxial  tension  when  m = 1.  However,  for  m = 5 and  10,  it 
underestimates  small  equibiaxial  failure  probabilities  by  factors  of  approx- 
imately 2 and  3,  respectively.  Note  that  eq.  (9),  can  be  approximated 
surprisingly  well  by  the  very  simple  equation 


in  Ps  (o,  a) 
in  Pg  (a,  o) 


£2  (m) 


0.45 


(10) 


From  Fig.  1,  it  is  clear  that  the  BF  approximation  is  conservative 
for  m < 1.  However,  values  of  m in  the  range  of  0 < m < l are  abnormal. 
With  the  use  of  eq.  (2),  it  is  readily  shown  that  the  slope  of  the  P^(o)  curve 
is  zero  at  a = 0 when  m>  1,  in  agreement  with  observation.  For  m = 1,  the 
slope  is  finite;  for  m < 1,  it  is  infinite;  such  behavior  is  rarely,  if  ever, 
encountered  in  nature.  Thus,  we  conclude  that,  for  practical  purposes, 
eq.  (1)  can  be  considered  unconservative. 

We  now  consider  the  way  the  size  of  the  error  in  eq.  (1)  varies  with 
stress  ratio,  again  focusing  attention  on  the  biaxial  case.  It  was  pointed  out 
earlier  [8]  that,  for  materials  that  obey  eq.  (2),  the  Weibull  theory  and  that 
of  Batdorf  and  Crose  [ 10]  give  the  same  results.  In  the  general  case  of 
triaxial  tension  characterized  by  principal  stresses  ct^,  a 2 and  a^,  where 
crj  > ct2-  ct3  - °*  bot;h  theories  imply  that  the  probability  of  fracture  takes 
the  form 
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in  Ps (ct,(t) /in  Ps  (<x, ol 


Pj.  = 1 - exp 


2 2 
cos  p sin  0 


+ 02  sin 


0 sin  0 


m 


+ o,  cos  0)  sin  0 d 0 


(ID 


where  0 and  0 are  the  polar  and  azimuthal  angles,  respectively.  As  noted 
by  Shur  [ll],  the  integral  in  eq.  (11)  yields 


sin  0 d6  = 4tt  o j 


(ml  )2 
(2m  + 1)J 


m 

y (2i)l  (2j)I  (2k)!  j k 
i,j,k  (il  j!  k')2 


(12) 


where  i,  j,  k are  any  positive  integers  that  satisfy  the  inequali' ies 
o < i,  j,  k < m and  i + j + k = m,  and  where 


s = (13a) 

t = ct-j/Oj  . (13b) 


For  the  biaxial  case,  o^  = t = o,  and  the  probability  of  survival  becomes 


where 


P 


s 


4rrk 
2m  + 1 


F (m,  s) 


(ml ) Y''  (2i)  ! [ 2(m  - i) ] l 

{2m)l  to  iiMm-i)!]2 


(14) 


(15) 
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Thus,  in  the  general  biaxial  case, 


in  Ps  (gl,  g2)  ~ F (m>  s) 
in  Pg  (oj , o)  " F (m,  o) 


F (m,  s) 


where  m is  any  positive  integer,  and  o < s = — < 1. 

al 


The  BF  approximation  gives,  for  the  biaxial  case, 


in  Pg  (cj,  o2)  in  Pg  (aj,  s) 

in  Pg  (aj,  o)  = in  Ps  (oj,  o)  = 1 + s 


(16) 


(17) 


Curves  that  represent  eqs.  (16)  and  (17)  are  shown  in  Fig.  2.  It  is 
evident  that,  for  values  of  rn  of  practical  interest,  the  BF  approximation 
consistently  underestimates  the  fracture  contribution  of  the  second  stress, 
and  that  the  error  is  largest,  in  absolute  terms,  at  least,  for  the  equibi- 
axial  case,  s = 1 . 

Another  way  to  compare  theory  with  the  BF  approximation  is  to  take 
the  ratios  of  eqs.  (16)  and  (17).  Thus, 


in  (Ps^Th  _ F (m,  s) 


*n  ^^Vbf 


1 + s 


m 


(18) 


This  result  is  shown  in  Fig.  3.  It  is  evident  that  the  ratio  increases  mono- 
tonically  from  s = o to  s=  1.  Moreover,  in  the  range  2 <m  < 10,  the  ratio 
is  reasonably  constant  near  s = 1. 

The  first  of  these  observations  suggests  a simple  method  for  obtaining 
an  upper  bound  to  the  failure  probability,  or  a lower  bound  to  the  survival 
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probability  of  a biaxially  stressed  structure;  viz.  , it  implies  that,  in  every 
element  of  the  structure, 


*n  (Ps}Th  . F (m,  s)  ^ F (m,  1)  ^ 0.45 

£n  ^Ps*BF  1 + sm  2 


Since  eq.  (19)  holds  for  every  element  individually,  it  must  hold  for  the 
structure  as  a whole. 

Thus , we  conclude  that 

- in(Ps)BFS  -^(Ps)Th<-m°-45i„(Ps)BF  . (20) 


For  small  values  of  P^.,  this  implies  that 

^Pf^BF  ~ ^Pf^Th  m ^Pf^BF 


(21) 


Because  the  curves  of  Fig.  3 are  fairly  flat  near  s = 1,  in  the  case  of 
structures  whose  loading  is  nearly  equibiaxial,  the  second  inequality  in 
eqs.  (20)  and  (21)  approaches  equality.  In  the  case  of  circular  disks  with 
uniform  lateral  or  concentric  ring  loading,  for  instance,  the  (small)  prob- 
abilities of  failure  should,  according  to  the  theories  discussed  here,  be 
0 45 

quite  close  to  m ’ times  the  result  given  by  the  BF  approximation,  which 
yields  a simple  closed-form  solution. 


As  mentioned  earlier,  the  preceding  analysis  is  based  on  the  assump- 
tion that  the  uniaxial  tensile  fracture  statistics  of  the  material  under  con- 
sideration can  be  described  with  adequate  accuracy  by  Weibull's  two- 
parameter  equation.  When  the  material  is  better  represented  by  his 
three-parameter  equation,  the  analysis  becomes  more  complicated  and 
is  beyond  the  scope  of  the  present  paper.  However,  the  general  observa- 
tion that  eq.  (1)  is  unconservative  still  applies,  and  it  is  to  be  expected 
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that  the  size  of  the  error  might  be  substantial,  particularly  when  the  Weibull 
parameter  m is  large.  Finally,  it  should  be  pointed  out  that,  for  materials 
that  obey  eq.  (2),  an  accurate  Weibull- type  treatment  of  structural  problems 
that  involve  biaxial  or  triaxial  stresses  is  not  difficult  when  a computing 
machine  is  used  in  conjunction  with  finite  elements  provided  that  closed- 
form  solutions  to  the  integral  in  eq.  (11)  are  used.  These  solutions  are 
given  in  eqs.  (12),  (14),  and  (15). 
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LABORATORY  OPERATIONS 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  arc: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromi  gnetics:  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  p iysics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Physics  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airgiow,  xudgneluspheiic  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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